Baryons and weak lensing power spectra 
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Abstract 

We provide a simple toy model to elucidate the effects of baryonic cooling on predic- 
tions of the convergence power spectrum in weak gravitational lensing. Our model 
suggests that the effects of baryonic cooling on dark matter halo profiles can al- 
ter the dark matter power spectrum at k ~ 10/iMpc -1 and lead to percent level 
changes in the convergence angular power spectrum at wavenumbers of several thou- 
sand. These effects would be measurable with planned weak lensing surveys and will 
impact constraints on cosmological parameters. 
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1 Introduction 



The study of modern cosmology has been tremendously advanced by probes 
for which detailed comparison of theory and observation is possible. The pro- 
totypical example is the anisotropies in the cosmic microwave background 
(CMB) for which t heoretical calculations have reached a highly ref ined level 
f Seliak et al. . 2003j ) and measurements a high degree of accuracy (e.g. Bennett et al 
2003|). Weak gravitational lensing, the small shearing of the shapes of back- 
ground galaxies by gravitational lensing from large-scale structure along the 
line-of-sight, is another. With detections of the effect now becoming possibl e 
and a relatively clean theory ( Mellier . 19991 : Bartelmann fc Schneider! . l200lh . 
gravitational lensing is beginning to become a critical cosmological probe, al- 
ready pr oviding constraints on the mass density and the fluctuation amplit ude 
(see e.g. Ivan Waerbeke fc Mellierl . 120031: iHoekstra. Yee. fc Gladder! 120021 for 
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the current status) and touted for its potential to constrain other interesting 
cosmological parameters. 

One of the reasons often quoted for the cleanness of weak gravitational lens- 
ing is that the deflection of light involves only gravity, and this circumvents 
many of the intractable problems associated with complex baryonic physics 
(e.g. cooling, feedback, star formation etc.). Clearly this statement is an ideal- 
ization. The ability of baryons to cool and collapse into dense structures alters 
the gravitational field around them and will thus alter the predictions of grav- 
itational lensing. The dense cores of galaxies can exceed the critical density 
necessary for the production of multiple images, strongly modifying the tail 
of the lensing distribution. The cooling of baryons also causes a concentration 
of the halos in which they live, altering the matter power spectrum on small 
scales. 

We wish to investigate this latter effect in this paper, using an analytic toy- 
model to gain a feel for the expected amplitude. This will allow us to estimate 
the "theory uncertainty" associated with the predictions, and provide guidance 
on experimental design. As we shall see, precise measurements of the small- 
scale lensing power spectrum will be sensitive to the details of baryonic physics, 
providing another interesting probe of these complex phenomena. 



2 Relevant scales 



We begin by discussing the relevant length and angular scales before turning 
to our toy model for including baryonic cooling. We shall work throughout 
in terms of the power spectrum of the convergence, k. This, or equivalent 
measures, has received much of the attention from the community to date. 
The convergence is a weighted integral of the mass density along the line-of- 
sight: 
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where S is the overdensity, a is the scale-factor, x is the comoving distance 
and g(x) is the lensing weight 
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for sources with distribution p(Xs) normalized to J dp = 1 
and Born approximations we can write the (di mensionless) 
k as a n integral over the mass power spectrum (jMellierl ^ 
2001) 



(2) 

Using the Limber 
ower spectrum of 
iBartelmann fc Schneider . 



9vr 

Ti 



O 2 ff 4 

"mat^O 



X' dx 



A 2 m (k = £/ X ,a) 



(3) 



2 



where g(x') is defined in Equation 2, A^(fe) = k 3 P(k) / {2ir 2 ) is the con- 
tribution to the mass variance per logarithmic interval in wavenumber and 
A 2 K (£) = £ 2 Ct/(2iT) is the contribution to 

K rms P er logarithmic interval in an- 
gular wavenumber I. 

Our ability to measure A 2 is constrained on both large and small scales. On 
large scales the finite amount of sky we can cover leads to a "sample variance" 
error. For a Gaussian k field, a very good approximation on scales above a 
degree, this leads to a fractional error 
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where / s k y is the fraction of sky surveyed. On small scales we are limited by 
the finite number of galaxies with which we can measure the shear. If n ga \ is 
the number density of galaxies of intrinsic ellipticity 7,-ms, this shot-n oise term 
increases the error to (|Mellierl . Il99flt iBartelmann fc Schneider! l200lh 



5 °± = / 2 (l I 7 ™ ^ (5) 

C t V (2£+ l)/sky V %alCV 

Ambitious experiments now being planned hope to measure A 2 over the range 
10 2 < £ < 10 4 using many tens of galaxies per square arcminute over much 
of the sky 2 . For an all-sky experiment with n g3i \ — * oo we could measure the 
power in a 10% band around £ = 10 3 to 0.3%. This motivates us to consider 
what the theoretical uncertainty on our predictions is. 



3 A model 



Previ ous work (e.g. lJain. Seliak fc White] . l2000HVale fc WhitelEooilWhite fc ValeL 



2004) has studied the level of accuracy of assumptions in the Limber and Born 
approximations and numerical artifacts arising from the calculational proce- 
dure under the assumption that only gravitational physics operated. Here we 
want to consider the effects of baryonic physics, specifically cooling, on the 
power spectrum. To set the stage we show in Figure 1 how a (redshift inde- 
pendent) change in the matter power spectrum translates into a change in 
the angular power spectrum of k through Equation 3. We have artificially 
increased the small scale power by multiplying A 2 nsBS by 1 + (k/kj) 3 for a 
range of kf. Note that setting kf = 30/iMpc _1 corresponds to a 4% increase 
in power at k ~ 10/iMpc _1 , while setting kf = 10/iMpc -1 is a factor of 2 



2 see e.g. jht tp: / / pan-starrs. org/| , |http://www.c tio. noao.edu / telescopes / dec. html 
http://www.lsst.org and http://snap.lbl.gov 
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Fig. 1. The effect of (artificially) increasing the small scale mass power spectrum 
on the angular spectrum of n. The top panel shows the angular power spectra for 
sources at z s = 1 (dashed) and 2 (dotted) when A^ ass (A;) is multiplied by l + (k/kf) 3 
for kf = 00, 30, 10 and 3/iMpc -1 (in order of increasing small scale power). The 
lower panel shows the power ratio, compared to kf = 00. 



increase on the same scale. The comparison of the shift in the signal to the 
estimated uncertainty is given in Figure 2 for sources at z s — 1. 



At what scale do we expect baryonic physics to be important? The dominant 
effect of hot baryons is to isotropize the inner regions of halos, due to the finite 
(isotropic) pressure they experience. This has a negligible effect on the power 
spectrum. Thus we are more interested here in the baryons which can cool. 
We shall assume that the dominant effect of cooling is to cause the baryons 
to contract, potentially dragging some of the dark matter with them through 
their gravitational coupling. Since substructures make up a small fraction of 
the total mass of a halo we shall consider only the baryons cooling to the 
center of any given halo, which should dominate the effect. 



To gain a feel for the size of the effect of cooling on the mass profiles of halos 
we make use of the adiabatic cooling model. We assume "sphericized" halos 
containing dark matter and hot and cold baryons. We ignore the pressure of 
the hot baryons, treating them in the same m anner as the (cold) dark matter . 
The initial halo is modeled as an NFW profile (jNavarro. Frank fc Whitel Fflflfi) 
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Fig. 2. The fractional error in Cg, vs. multipole I for sources at z s = 1. The solid 
line shows the shift in Ct for kt = 30 /iMpc -1 from Figure 1. The results of section 
4 suggest this is a good ballpark estimate for the effect of baryons. The dashed 
lines show the cosmic variance limit for Gaussian fluctuations in a band of width 
A£ = 0.1 £ for / s k y = 0.1 (upper) and 1 (lower). The dotted lines include shot noise 
for galaxy number densities n ga i = 100, 50 and 25 galaxies per square arcminute 
(bottom to top) assuming the fiducial model, kf — > oo. 

with 

M(< r) oc log(l +x) - — — (6) 

1 ~~ r - X 

where x = crj r v - a is a scaled radius for a h alo of concentrat ion c. The baryons 
are assumed to follow a Hernquist profile ( Hernquist . 1990t 



x 2 



M(< r) cc — (7) 
(1 + x) z 

with x = r/a. Half of the mass is contained within r ~ 2a in this model. 

For a massive halo, say a rich cluster of galaxies, the concentration is near 5 
and the virial radius is ~ 1.5 h^ 1 Mpc while the cooling radius is O(100 kpc) so 
we take the initial scale = 0.05 — 0.1 r vir . Cooling usually allows the baryons 
to contract by a factor of 10 in linear scale, so we take af = 0.1 a^. Our final 
results are quite insensitive to the choice of within this range, and even to 
the choice of the other parameters to a large extent. 

Adiabatic contraction conserves both the mass and the angular momentum of 
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Fig. 3. The effect of adiabatic cooling on NFW halos. The change in the enclosed 
mass is plotted versus radius. The solid lines are for / coo i = 5% while the dashed 
lines are for / coo i = 10%. The lower curves in each case have c = 5, en = 0.1 r v ; r and 
a/ = 0.1 aj while the upper curves have c = 10 and a, = 0.3r v ; r . 

the co ld baryons and dark matter particles ( Blumenthal et al. , 198fit Kochanek fc White! 
2001). Thus given a final radius we can solve for the initial radius through 



M f {r f ) = M i {r i ) + M cf {r f ) - M Ci (r,) 
r f M f (r f )=r i M i (r i ) 



(9) 



where Mj and Mf are the initial and final total mass profiles and M C j and 
M c f are the initial and final distributions of the material which can cool. We 
show the effects of this in Figure 3. The cooling and contraction of the baryons 
draws material into the center, enhancing the central density. The effect is not 
very sensitive to the initial halo concentration and scales almost linearly with 
the cooled mass fraction. 



4 Effect on mass power spectrum 



To see what the effect o f these changes are on t he matter power spectrum we 



use the halo model (see ICoorav fc Shethl . I2002L for a review). While this will 



not serve for high precision cosmology, it will illustrate the typical scale of the 
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effects. The fundamental assumption in this model is that all of the mass in 
the universe resides in virialized halos whose clustering properties are those of 
peaks of the density field and whose profiles are drawn from a universal one 
parameter family depending only on halo mass. The power spectrum is then 
made up of two contributions 

P{k) = P 2 - halo (A;) + P^^ik) (10) 

one from pairs of mass elements which lie in different halos and one where they 
both lie in the same halo. The first term dominates on large scales, while the 
second dominates on the smaller scales of interest to us. Within this model we 
can compute the mass weighted average of any quantity by integrating over 
the multiplicity function 

where v is the peak height. For f (v) we use the fit to N-body simulations 
quoted in (|Sheth fc Tormenl . ll999h . 



The power spectrum from the 2-halo term is that due to a system of (smooth) 
halos of profile y(k) laid down with inter-halo correlations assumed to be a 
biased sampling of Pn n (k). This dominates on large scales. Since the real space 
convolution is simply a Fourier space multiplication this contribution is 



2 

->2— halo / 



P 2 - hal °(k) = P lin (k) 



f(u)du b(u)y(k;M) 



:i2) 



where b{v) is the (linea r) bias of a halo of mass M{y) which we take from 
( Sheth fc Tormenl . E999). We assume that the halos all have spherical profiles 
depending only on the mass. We neglect any substructure or halos-within- 
halos. We take the NFW form with Fourier transform 



p(k) = 47rpo^ [ cos2{Ci([l + c]z) - Ci(z)} + 

sin cz 



sin,z{Si([l + c]z) -Si(z)} 



(13) 



[l + c)z_ 

where z = kr s and Ci and Si are the cosine and sine integrals. 

On small scales we are dominated by pairs lying within a single halo 



where p is the mean matter density. This term gives the major contribution 
to the angular power spectrum of k on the scales of interest. 
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Fig. 4. The effect of cooling baryons on the mass power spectrum. The 
solid line is the unmodified power spectrum, the dotted line includes the 
correction for e = 0.05tanh(10 14 h~ 1 M Q / M V - 1T ) and tq = 0.05r v i r , the 
dashed line for e = 0.1 tanh(10 14 h~ 1 M< 3 /M v i r ) and the long-dashed line for 
e = 0.2tanh(10 14 /i _1 M /M v i r ) to illustrate the sensitivity to our assumptions. 
The lower panel shows the ratio to the unperturbed result. By z ~ 1 the increase is 
shifted to the right by a factor ~ 1.2. 

Within this formalism it is easy to see how we incorporate the effects of cooling 
baryons. The 1-halo term, which dominates the scales of interest, is simply 
modified by y(k) — > y(k) + 5y{k) where 



for changes in the cumulative mass, AM(< r), which vanish at r = and 
infinity. Here ji(x) is the spherical Bessel function of order 1. 

In keeping with our attempts to construct an analytic model of the effect we 
fit AM (< r) with a simple form 




(15) 




(16) 



which is a reasonable fit to the curves of Figure 3. Then 5y{k) is 



5y(k)=e-(kr ) [l - (1 + kr ) e~ kro 

= e(n/4:)( y kr ) 3 + --- kr < 1 



(17) 
(18) 



allowing us to compute 5P(k) as a function of e and r . It is also possible to 
find an analytic transform for (r/r ) 3 /(l + (r/r ) 4 ), which has a flatter core 
as r — > 0. The large-r and low-A; behavior is the same and for the /c-range of 
interest the results are almost indistinguishable. 

Fitting the massive halo differences shown in Figure 3 we find e ~ 0.15 and 
r — 0.05 r vir . We expect the cooled fraction to be very small in high mass 
halos and to increase with mass, until it reaches a limit set by fibar/^mat- 
Taking e ~ 3 / coo i we model this as e = (0.03 - 0.1) tanh(10 14 h^M^/M^). 
We also expect ro to become a larger fraction of r v - 1T for lower mass halos - 
the trend is weak however, so we hold ro = 0.05 r v ; r for simplicity. The precise 
details of the low mass behavior are unimportant for us, as it is the most 
massive halos that dominate the k ~ 10 h Mpc -1 regime which will be our 
focus. Also for kr <C 1 there is a degeneracy between the effects of a change 
in r and a change in e. 

Using the e and r above gives the power spectrum shown in Figure 4. In this 
model there is a few percent increase in power at k ~ 10 h Mpc -1 and a factor 
of two increase at k ~ 20/iMpc" 1 . By z ~ 1 the curves have shifted to the 
right by a factor 1.2. The ratio is reasonably well fit by the 1 + (k/kf) 3 form 
used earlier with kf ~ 30 /iMpc -1 . Thus Figure 1 suggests this would lead to 
a few percent increase in the power at £ ~ 3000 — 5000, which is potentially 
measurable in upcoming experiments (see Figure 2). There is at least a factor 
of two uncertainty in our estimate of this scale, but it shows that we expect 
percent level modifications to the matter power spectrum at k ~ 10/iMpc -1 
and the convergence power spectrum at £ ~ 3000 — 5000 from baryonic cooling. 

Currently we are unable to reliably compute the effects of cooling, star for- 
mation and feedback on the small-scale mass distribution. This uncertainty 
will need to be folded into our requirements on theoretical calculations of the 
spectrum, or mitigated in some way. One possible route to reducing the sensi- 
tivity to the uncertain high-/c power is to make use of multiple source redshifts. 
Within the Limber approximation of Equation 3 only modes transverse to the 
line-of-sight contribute. This means that at fixed wavenumber, £, higher red- 
shift corresponds to larger wavelength. By combining signals from a number 
of different source redshifts the sensitivity to low- 2 structure can be reduced 
("nulling tomography"), and with it the sensitivity to high- A; at any given £. 
The price is a decrease in signal-to-noise. Since the number of modes on the 
sky depends on £, such a strategy would allow more accurate measurement 
of the large-scale power spectrum, uncontaminated by the higher k modes, if 
rtg a i was sufficiently large. Alternatively we could attempt to model the effect 
of baryon cooling. The total theoretical uncertainty is then reduced to the 
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uncertainty in our modeling. The treatment above is fairly crude but a more 
sophisticated modeling, informed by numerical simulations, may be possible. 



5 Conclusions 

We have considered the effect of cooling of baryons on the small scale mass 
power spectrum. With a simple model for adiabatic compression we can use 
the halo model to estimate the size of this effect on the 3D mass power spec- 
trum, and hence on the small angular scale convergence power spectrum for 
weak lensing. Our model contains numerous simplifying assumptions and pro- 
vides at best a "factor of 2" level calculation, but we find that baryons can 
alter the mass profile non-negligibly on scales of 100 /i _1 kpc which leads to 
modifications of the angular power spectrum of k at i ~ 3000 — 5000 at the 
several percent level. As indicated in Figure 2, these modifications could be 
probed by future weak lensing surveys and this would give us further insight 
into the inner parts of dark matter halos. It would also serve as an extra 
source of theory uncertainty in constraining cosmological parameters unless 
observational strategies to minimize it were implemented. A refinement of this 
calculation, using modern numerical simulations which model the heating and 
cooling effects of baryons, is clearly in order. 

MW. thanks Gary Bernstein, Dragan Huterer, Eric Linder, Chris Vale and 
David Weinberg for useful conversations and comments on an earlier draft. 
This research was supported by the NSF and NASA. 
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